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Abstract 

A  frequency-domain  stability  criterion  is  pre¬ 
sented,  generalizing  the  well-known  circle  sta¬ 
bility  criterion  to  maltiloop  feedback  systems 
having  bounded  nonlinearity ,  parameter  varia¬ 
tions,  and/or  frequency -dependent  ignorance  o: 
component  dynamics.  Unlike  previous  general¬ 
izations,  the  theory  is  not  restricted  to  weakly- 
coupled,  diagonally  dominant  or  nearly  normal 
systems.  Potential  applications  include  the 
analysis  of  feedback  .system  integrity  and  multi¬ 
loop  feedback  system  stability  margins,, 

1 .  Introduction 

A  key  step  in  the  synthesis  of  robustly  stable 
feedback  systems  is  the  characterization  of  a 
set  of  feedback  laws  that  are  stabilizing  for 
every  element  of  the  set  of  possible  plant  dy¬ 
namics.  This  type  of  information  is  precisely 
what  is»  provided  for  single-loop  feedoack  sys¬ 
tems  by  such  input -output  stability  criteria  as 
the  Nyquist,  Popov,  and  circle  theorems.  In¬ 
deed,  the  practical  merit  of  classical  feedback 
design  procedures  involving  Nyquist  loci.  Bode 
plots,  and  Nichols  charts  is  in  a  large  measure 
directly  attributable  to  the  fact  that  these  design 
procedures  provide  the  designer  with  an  easily 
in^erpretable  charcterization  of  such  sets  of 
robustly  stable  feedback  laws.  Available  multi- 
variable  input -output  stability  criteria  such  as 
Rosenbrock's  multivariable  Nyquist  theorem  f  1  ] 
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and  Zanies’  conic -relation  and  positivity  stability 
theorems  [2]  have  led  to  useful  characterizations 
of  sets  of  robustly  stable  feedback  laws  for  only 
a  limited  class  of  problems,  viz.  interconnec¬ 
tions  of  dissipative  systems  [3],  weakly  coupled 
interconnections  of  systems  [4]  (including  so- 
called  ’’diagonally  dominant”  systems  [5]  -  [7]) 
and  ’nearly  normal"  systems  [8]  (which  can  be 
viewed  as  vector-space  isomorphisms  of  weakly 
coupled  systems).  It  has  been  argued  convinc¬ 
ingly  by  Rosenbrock  and  Cook  [9]  that  an  es¬ 
pecially  useful  feedback  design  tool  would  be  a 
more  general  multiloop  frequency-domain  sta¬ 
bility  criterion  that  includes  diagonal  dominance 
and  normality  results  as  special  cases. 

The  main  result  of  the  present  paper  is  a 
stability  result  that  may  serve  this  purpose; 
Theorem  l  is  a  multiloop  generalization  of  the 
circle  stability  criterion  which  does  not  re¬ 
quire  diagonal -dominance,  weak-coupling,  norm¬ 
ality,  or  near  normality.  The  result  allows  the 
frequency-domain  testing  of  the  stability  of 
muitiloop  feedback  systems  with  time-varying 
nonlinearities,  unknown-but -bounded  parameter 
variations,  and  even  singular  perturbations. 

T  T 

The  following  notation  is  used;  A  and  x 
denote  respectively  the  tjranspo|e  of  the  matrix 
A  and  the  vector  x  ;  A  and  x  denote  the 
complex  conjugate  of  the  matrix  A1  and  the 


complex  conjugate  of  the  matrix  A  and  the 
vector  xA  respectively;  the  determinant  of  a 
matrix  A  is  denoted  dcr(A);  fhe  Euclidian  norm 
of  a  vector  x  is  jx  =^x  x  ;  R^  denotes 
nonnegative  real  numbers;  the  functional  norm 
.  x  (<T  and  inner  product  <<^,  x?  >-  are  defined 
tor  functions  x:  ■*  Rn  as 


x  _  =  v<  X,  X  >_ 


where  for  any 


le  space 


J0  lt>  x,  (t)  dt. 

)  is  defined  as 

’  {*  =  R_  .<• 

for  each  ~tR  >  . 
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Given  any  matrix  A,  the  square -roots  of  the 
eigenvalues  of  A"  A  are  called  the  singular  values 
of  A  [10,  pp.  5-11],  For  any  non- zero  Matrix  A,  we 
use  the  notation  ama^A)to  denote  the  largest  sing¬ 
ular  value  of  A  and  ami^A)to  denote  the  smallest 
singular  value  of  A;  singular  values  are  always  non¬ 
negative  real  numbers  since  A*  A  is  always  posi¬ 
tive  semidefinite.  For  hermitian  A  (L  e. ,  A  =  A*), 
the  notation  ^^A)  and  Vnin^  t0  denote  respec¬ 
tively  the  greatest  and  least  eigenvalues  of  A - 

hermitian  matrices  have  only  real  eigenvalues, 
so  ordering  of  eigenvalues  is  always  possible  via 
the  usual  ordering  of  real  numbers. 

An  operator  is  a  mapping  of  functions  into 
functions;  for  example,  a  dynamical  system  map¬ 
ping  inputs  x  ^  L-2e  (R+*  R  )  into  outputs  y^L^e 
'R+,  Rm>  defines  an  operator.  All  operators  con¬ 
sidered  in  this  paper  are  implicitly  assumed  to  be 
mappings  of  L^e  (R+,  Rnl )  into  L^e  (R+,  R°2) 
for  some  positive  integers  ana  n ?  #  an  opera¬ 
tor  His  said  to  be  nonantic ipative  if 


Dynamical  LTI  interconnection 


(  H  xt)  (t0 


<  5  *2  >  (to 


for  any  tQ  and  any  pair  of  functions  and  x, 
having  the  property  that  for  ail  t  <  t  , 


i.e,,  a  non  anticipative  operator  H  is  one  having 
the  property  that  its  instantaneous  output  Hx  at 
any  time  t  is  independent  of  the  values  assumed 
by  the  input  x(t)  at  future  times  t  >  t  •  We 
say  that  an  operator  H  mapping  signals 
x  £  L^e  {R^,  R"  )  into  signals  Hx  €  L^e  »R+,  R™) 
is  L2e  -  stable1  if  there  exists  a  constant  k 


X(s)  =  -  G(s)(y(s)  +  V(s))  +  U(s)  (2.3) 


Y(s)  £ 


;  x(s)  r 


Y  (s) 
.  ntm 


:U(s)S 
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4  INTERCONNECTION 


I  :i 

Ji  il 


such  that  for  all  x  €  L, 


R  I  and  tcR 


II.  Problem  Formulation 

Our  results  concern  the  input-output  stability 
of  systems  consisting  of  a  dynamical  LTI  nega¬ 
tive-feedback  interconnection  of  m  memoryless 
time -varying  nonlinear  components  and  n  dyn¬ 
amical  LTI  components.  The  system’9  equa¬ 
tions  thus  take  the  following  form  (see  Fig.  L): 
Memoryless  nonlinear  components 

y.  <t)  =  h  {x  (t)  ,  t  )  i 

:  (2.D 

y  (t)  =  h  (x  (t)  ,  t  )  > 

m  mm 

Dynamical  LTI  components 

Y  .<S>?  I  <*>  X  „<*>  I 

m  1  ;  m'1  mfn  (2.2) 

Y  (s)  =  H  (s)  X  (s)  S 

_  ...m.rn.  t  _ _ _ 

1  A  form  oi  global  input-output  stability,  our 
definition  of  L2e  *  stable  is  equivalent  to 
’  L  stable  1  [11]  and  to’finite  gain  stable1 
{with  respect  to  the  I.,  norm)  [12]  when  H 
is  a  nonanticipative  operator. 


Fig.  1 _ The  System  _ 

The  indogenous  variables  y^(t)  and  x.(t)  are 
the  system  ’outputs’  and  the  exogenous  variables 
u^(t)  and  v.(t)  are  the  system  ’inputs'.  Each  of 

the  ’components’  (h^(»  ) . hm,  Hm+  p  .  . . , 

Hmfn)  may  itself  be  a  M1MO  system  in  general, 
though  our  results  seem  to  be  most  easily  used 
and  interpreted  when  the  components  are  SISO. 

Our  stability  results  do  not  require  that  we 
have  available  an  exact  mathematical  description 
of  the  components.  For  each  of  the  nonlinear 
nondynamical  elements,  we  assume  only  that 
matrices  C.  .  R.  t  and  5.  can  be  found  such  that 
R.  R.  and  ^.^S^are  positive  definite  and  suchtiiat 

))  -  s 

for  some  «  >  0,  all  x^(t),  and  all  t 
( i  =  1  f  .  .  .  ,  m), 

(2.  5i 

For  each  of  the  n  dynamical  LTI  components 
H4(s>  (  i  =  m*I . m-n)  we  assume  onlv 


MILL'  'H'wq 


tha*  IL(s)  has  a  proper  rational  transfer  function 
matrix  and  that  proper  rational  transfer  function 
matrices  C^(s),  R^(s),  and  Sj(s)  can  be  found  such 
that  Rj  (-jin)  R^jjj  )  and  Sj^(jju  )  S^ju)  are  positive 
definite  and  have  no  poles  on  the  3  =  ji;  axis  and 
such  that  H^s)  -  C.(s)  has  no  pole3  in{  s  |  Re(s)^oJ 
and 


1  S.(j»  ^(j* )  X.tjai  )  -  C.(j*  >  X.(j»  |ji!2 
*  ;  R.fjjox.o*)  ;| 2  - .  .|x.(jji)!|2 


(2.6) 


for  some  «  >  0,  all  X^(jji)  and  all  ui 


(i  =  m+1 . m+n). 


For  notational  convenience,  we  define  the  follow* 


ing  block -diagonal  matrices 

C(s)  =  diag^C ^ , , . . 

’Cm’  Cm+l(S)-  ” 

. .  C  (sA 

m+n  ) 

(2.7) 

R(s)  =  diag^Rj , .  . . 

••Rm+H 

(2.8) 

S' s )  =  diag^Sj , . . . 

,  S  ,  S  (s),.. 

m  m+1 

•'Sm+n(S)) 

(2.9) 

• — h  (,*5  ) 


Comments: 


The  conditions  (2,5)  -  (2.6)  can  be  interpre¬ 
ted  as  saying  that  we  are  assuming  knowledge 
about  each  of  the  components  is  limited  to  an  ap¬ 
proximate  LIT  model  (viz.  C.J  and  bour.dsfde- 
termined  by  (R.  ,  S.))  on  the  coarseness  of  the 
approximation. 


For  the  case  of  SISO  components,  the  condi¬ 
tions  (2.5)  and  (2,6)  can  be  replaced  respectively 
by  the  simpler  conditions 


|h.(x.(t),t)  -  Cix.(t)|2 


for  some  *  >  0  ,  ail  x.(t)  0,  and  alL  t 

|  H.(jn )  -  c.UxlfsIr.ljjof  -  . 


lor  some  «  >0  and  ail  u 

where  for  all  i  =  1,..,,  n+m 

c, 
t 

r. 

1 


R.  S. 

I  L 


-l 


(2.5«) 

(2.  6’) 


(2.  10) 
l2.  11) 


(b)  Nyquist  locus  of  LTl  component 
satisfying  (2.6*) 

Fig-  2  51SQ  Components _ 

II I,  Main  Result 


Our  main  result  is  now  stated. 

Theorem  I  (Multiloop  Circle  Stability 
Criterion); 

Suppose  that  the  system  (2*  l)-(2.4)  is 
L2e-*table  for  the  case  when 

h.(x,t)  *  C.x  (i  =  1 . m)  \ 

(3.1) 

\\{s)  E  C^s)  (Urrw-1,...,  mm),  J 


i) 


Then,  a  sufficient  condition  for  the  system 
(2,1)  -  (2.4)  to  be  L2e  -stable  for  every  col¬ 
lection  of  h.(  •  .  t)  (i  =  1 ,  . .  „m)  and  H;(s) 

(i  =  m+1, . . . ,  n)  satisfying  (2.  5)  and  (2.  6) 
respectively  is  that  any  one  of  the  following 
conditions  hold  for  all  real  w 


\inj^I+CfjJU  IGLjju  )jr  ;  S(j u )(l+C(  jau  )G(j  *)) 


-G  (-joj)R  f-jju  G(jJ>)]  *  0 


(3.  -a) 


ii) 


\nii/(CW*  J+G'Vj*  S?-jx  lS(jx  )(c(j*>*GJ(j*j 
-R*  (-j-»  (  R(jt  ^  0 


These  31SO  conditions  are  readily  interpreted 
graphically  as  shown  in  Fig.  2. 


3 


(3.  2bi 


lU>  *  min  («j.  )  (  c (>  )  +  G'  l(j*  ))  R  1(jju  )j  a  1 

(3.2c) 

ivl  e'max(R(jA' )  G<**  *  (r  +  CO^O^S'W*) *  1 

(3.  2d) 


Condition  (3.2a)  is  implied  by  (3.  2b)-(3,  2d)  and, 
when  the  inverses  G'1,  R“ are  defined, 
conditions  (3.  2a) -(3.  2d)  are  equivalent. 


+ 


where 


l 

G(jju) 


r) s  r.s;‘ 


c.  =  ci 


<* 


or.  equivalently  (assuming  >  0), 


PROOF :  See  Aopendix 
IV.  Discussion 


There  are  essentially  two  main  conditions 
which  must  be  satisfied  to  conclude  stability  from 
Theorem  Is  (i)  The  system  must  be  stable  when 
the  uncertain  nonlinear  components  h^(«,t)  and 
1-TI  components  H^(s)  are  all  replaced  by  the  re¬ 
spective  LTI  approximations  and  C^(s);  and, 
(iij  the  frequency-domain  condition  (3.2)  must  be 
satisfied.  The  former  condition  can  be  verified 
a  variety  o£  svays:  for  example,  one  may  check 
that  the  roots  of  the  characteristic  equation 


dei  (i  +  C(s)  G(s)j  =  0 


(4.  1) 

all  have  negative  real  parts’";  alternatively,  one 
may  apply  the  multivariable  Nyquist  criterion, 
checking  that  the  polar  plot  of  the  locus  of 
de-  (C(j»)  G(jn ))  encircles  the  point  -1  +  jO 
exactly  once  counterclockwise  for  each  unstable 
open-loop  polecfc_(s)  G(s)  (multiplicities  counted) 
[1]  ,  [21],  The  latter  condition  (3.2)  requires 
that  one  plot  the  variable  a  (• )  or  a  (.  ) 
verses  i;  and  verify  that  th e  appr opr iate^ne qual¬ 
ity  holds  for  all  j.  . 


In  the  special  case  in  which  there  is  a  single 
scalar  nonlinearity  h^x^.t)  (so  that  m  =  l  and 
n  s  01#  both  of  the  conditions  of  Theorem  i  can  be 
verified  by  inspection  of  the  polar  plot  of  Gtj«) 
vs  4,  .  Stability  for  the  special  case  h^x^.t)  = 
is  assured  by  the  Nyquist  stability  criterion  if 
and  only  if  G(jx)  encircles  -1  +  jO  once  for  each 
unstable  pole  of  G(s)  as  jj  increases  from  -  »  to 
-  -  .  Condition  (3.  2)  becomes 


2  I:  C(7i  Gfs)  ha^  any  ‘decoupling  zeroes1  Vi.  e.  , 
uncontrollable  or  unobservable  poles),  then 
these  will  not  be  roots  of  *4.  1)  and  one  must 
check  separately  that  these  poles  have  negative 
reai  parts  -  ci.  (14] 


(4.3a) 


(4.  3b) 


(4.3c) 

These  conditions  on  G(j*)  are  precisely  the  con¬ 
ditions  of  the  well-known  circle  stability  cri¬ 
terion  (cf.  [15]  ).  It  is  this  which  motivates  us  to 
refer  to  Theorem  1  as  a  ’multiloop  circle  stability 

criterion'  - despite  the  fact  that  in  general  no 

circles  are  employed  in  verifying  its  conditions. 

One  can  interpret  the  uncertainty  bounds 
(R.f  S.)  as  specifications  for  the  gain  margins 
anci  pfiasc  margins  cf  the  system  (2.  1 )-(?..  4). 

If  rn»0,  if  H^s)  (la l . n)  and  if  the 

components  are  SISO,  then  under  the  conditions 
of  Theorem  1#  the  system  will  remain  stable 
despite  variations  in  the  individual  component 
gains  of  magnitudes  as  great  as  |r;(ju^  I  A 
|R,;  ( j  u)/S^  <  j  w)  |,  even  when  the  variations  occur 
3imuitac3ously  in  all  components.  So,  lor  ex¬ 
ample,  the  system  can  tolerate  simultaneous 
gain  variations  or  phase  variations  of  at  leas" 

T  UU>) 

°Ml  -  '•  "h  (4-4’ 

or 

^  r.ijui) 

3  inf  arcsin  |  — r-r  —  |  (4.5) 

M  4)  1  u.  (jx) 


i)  if  c 


I  G(jx ) 


1  1 
c , 


I  >  l 


2  C 1  \rl 

ii)  if  C,  -  T*  <0 


2  2 
C1  *  b 


G(jx  )  + 


'  1 


l  <  I 


...  2  1  2  1 
n.)  if  Cj  -  r,  =  0 


Re 


(GO"  )) 


b  *CI 


C1  +  1  rl 
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in  each  of  the  respective  component  input  chan¬ 
nel'  s  (i=l# . . . ,  n);  i.  e. ,  the  system  has  gain  mar¬ 
gins  of  at  least  ^d  phase  margins  of  at 

least  0^  at  the  inputs  to  the  respective  com¬ 
ponents  CL(s)  (i  =  l , . . .  , n).  The  quantity 


k  (c(j  U»)  +  G‘  \i  tu))  R  \j  -U)^ 


<4.  6) 


is  the  amount  by  which  the  uncertainty  bounding 
matrices  R^  can  be  simultaneously  increased 
without  violating  the  stability  conditions  of  The- 
crem  l  -—the  scalar  km.  can  he  viewed  as  a 
lower  bound  on  the  amount  by  which  the  system 
(2. l)-(2. 4}  exc eeds  the  stability  margin  speci¬ 
fications  (2.5)-{2.6*. 


In  general,  the  stability  conditions  ■ 
and  the  estimate  (4.  6)  of  excess  stability  mar¬ 
gin  k m  . of  Theorem  1  will  be  conservative. 

The  conservativeness  can  usuaily  be  reduced  by 
multiplying  equations  (2.  5) -(2.  6)  by  appropri¬ 
ately  chosen  positive  scalars  |  «.  |2  and 
|  arj(jo))  |  2  respectively  before  applying  Theorem 
1.  This  has  the  net  effect  of  substituting  'weight¬ 
ed*  uncertainty  bounding  matrices  (*i  R;  #  cr^S^) 
for  the  original  matrices  (R^  ,  ).  Further 

study  is  required  to  determine  the  extent  to  which 
it  will  be  practical  to  exploit  such  weighting  to 
reduce  the  conservativeness  of  Theorem  1. 


V.  Conclusions 


The  practical  importance  of  our  multiloop 
circle  theorem  is  that  it  provides  verifiable 
sufficient  conditions  for  the  stability  of  multi¬ 
loop  feedback  systems  using  only  crude  bounds 
on  system  parameters,  component  frequency 
responses,  and  nonlinearities.  Potential  appli¬ 
cations  include  the  testing  of  system  integrity 
in  the  presence  of  actuator  and  /  or  sensor 
failures  (cf.  [16]  )  and  the  characterization  of  the 
stability  margins  (e.  g.  ,  gain  6c  phase  margin)  of 
multiloop  feedback  designs  subject  to  simultaneous 
perturbations  in  gain  and  phase  in  the  feedback 
loops.  Thecrem  1  also  plays  a  key  role  in  bound¬ 
ing  the  response  of  systems  with  uncertain  dy¬ 
namics  f  22], 


Appendix 

In  this  appendix.  Theorem  1  is  proved  using 
the  sector  stability  criterion  of  [17],  [18],  [19]. 

We  begin  by  introducing  some  additional  termin- 
ology,  and  a  relevant  special  case  of  the  sector 
stability  criterion,  viz.  Theorem  Al.  We  then 
establish  via  several  lemmas  that  the  conditions 
of  Theorem  1  are  sufficient  to  ensure  that  the 
conditions  of  Theorem  Al  are  satisfied. 

DefinitionfL^g  -Cone<-,«,» );  strictly  inside,  outside): 


Given  any  three  operators  C,  R,  S,  we  define 

L2e  -  Cone  (C,  R,  S)  =  |  (x,  y)  J  F(x,  y,  t)  £  0 

for  aU  r  €  R  \  (Al) 
where  ' 

i)  =  il  S(y  -  Cx  II2  -  ||  Rx  ||2  (A 2) 

and  for  all  z,  z^ ,  z^ 

!!  *  :IT  =  V<  *.  z  >T  (A3) 

<  *r  z2  >r  “  /  ziT  dMA4> 

0 

Given  an  operator  H  mapping  signals  x  into 
signals  y  ,  we  say 

H  strictly  inside  -  Cone  (C,  R.  S)  (A 5) 

if  there  exists  an  «  >  0  such  that  for  every  pair 
(x,y)  satisfying  y  =  Hx 

Sis  (y-Cx)ilJ  S  ,  Rx  ii2  -  .(||  X  ||J+  It  y  ||2]  (A  6) 

for  all  T  €  R  . 

4. 

Given  an  operator  -G  mapping  signals  y  into 
signals  x,  we  say 

(-Gj1  outside  L.,^  -  Cone  (C.  R.  S)  (A7) 

if  for  every  pair  (x,y)  satisfying  x  =  -Gy  we  have 

il  S  <y-£*>  H*  a  ,|  Rx  !j^2  (A8> 

for  all  T  c  R  . 


Theor  em  A I  Let  p  be  a  positive  integer; 
let  H.  p)  be  operators  mapping  x  mt0 

y^;  let  H  be  the  operator  * 


II  s  <y-£x)  ilT2 

rT  2 

=  ^  il  s  (h(x(t).  t)  -  C  x(t)  ,j  it 

*  £  T||R  x(t)  II2  -  6  il  x(t)  II2  it 

*  ril-(t,ll2  -  e’(l|x(k)  ,|2 

+  iiy(t)  i|2  )  dt 

=  !|R*  ||t2  -  t  (  ||*  j|T2+  ■iyilT2  >  • 

(A  31) 

Conversely,  when  (A27)  holds,  then  taking  x(t) 
to  be  the  constant  function  x(t)  axQ  we  have 
that  for  some  6  >  0 

)|S  (h(*0  .  ^j)  -  C  ):i2 

=  -f-llA  (»*  -  Cx  )  L2 

*  — ~  (  il R*  il T2  -  S(H*liT2  +  !|y||T2^ 

S-i-(i|Rx||T2-4i|x|lT2) 

«  li  R*0  il  2  -  « II  x0  H2  •  <A3Z) 

cm 

Lemma  A3:  Let  H,  £,  J&,  be  nonanticipa- 
tive  L2e-*table  Unear-time-invariant  operators 
with  respective  rational  transfer  function  matri¬ 
ces  H(s),  C(s),  R(s),  S(a),  Suppose  that  R~  (s) 
exists  and  has  no  poles  in  Refs)  2  0* 

Then 


S  strictly  inside 
if  and  only  if 


L  —Cone  (C  ,  R,  S  ) 

2e  (A  33) 


li  s  (jo.)  (H(jK)  -  C(ji)))x  (jil)  il2  s  il R(jJJ  X (Kil 2 

-e  ii  X(jvu)  il2  (A34) 
for  aU  X(j^  all  w,  and  some  C  >0 


Proof:  Let  R~  denote  the  stable  nonantici- 
pative  LlTI  operator  having  transfer  £unction  ma¬ 
trix  R  (s).  Suppose  that  (A34)  holds  and  let 


f(Rx)(t)  ,  if  t  *  r 

zjt)  =  <  ~  (A 35) 

l  0  ,  if  t  >  t 

and  let  ZJ denote  the  Fourier  transform  of 
z^t).  Then,  for  all  y  =  H  x  we  have 

,S(y  -  Cx)1;2  =  S (Hx  -  Cx  )!lj 
=  ,1  S  (H  -  C  )  X 
(by  linearity) 

=  ii  S(H  -  C)R-‘rx  ' 2 
(since  R**  exists) 


»  il  S  (H  -  C)  R*1  xT  IL2 
(by  nonantic ipativea ess) 

*  jT il  (S(H-  C )  R1  *T»  (t)  :|2  dt 
(the  integral  exists  since  S,  H,  C,  and  iC^are  stdle) 

m 

*  f  ilS«(H(ja*-C(j<4>R‘l0u> 

||2  dm 

(by  Parseval’s  Theorem) 

*  £  I* 

dm 

(by  (A34)) 


=  f\ 'I'Jttlfdt 

ZS>  r  ,  , 

lies'  vc‘>i!  * 

(by  pars  aval*  s  Theorem) 


where 


*  f  ik(t>ii2<it 

'*•)  T 

-tjT  ii  (Rlv<‘>  ;i2  dt 

-  II  ’A'  -e  UsSC 

•  il  R*  IL2  -  6  ii  R"1  a*  iiT2 

(by  nonanticipativeness  of  rC*) 
=  il  rx||t2  -e  ;|  x  i^2 

‘  II R*  l2  -  c '  ( 11  x!|t2t||hx;|t2). 

(A36) 

* ' =  ~rh?  (A37) 


and 


Or 


=  sup 
x,  T 

(since 


(  ;l  8x  ik  \«r  • 

\  H*ilr  ) 

(A  38) 


H  is  stable) 


Conversely  suppose  (A 33)  holds.  Let 
and  ajj  be  arbitrary.  Then,  letting  x(t)  — ♦ 

e  and  t  •  f  we  have  trivially  from 
Parseval’s  Theorem  that 
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M  S(j^)  ( H(j-Ug)  -  C(j«tj))x0  II 

s  ilmj-t,)  X0  || 2  -  £  (  II  X0  if +li  H(jaij)  X0  II2  ) 

silR(J*b)  X0  ',i2  -  £  ||X0  i|2  .  (A39) 

CD 

Lemma  A4:  Let  G,  C,  K,  S  be  linear -time- 
invariant  operators  with  respective  proper  ration¬ 
al  transfer  functions  H{s),  C(s),  R(s),  S(s).  Sup¬ 
pose  that  S"  (s)  exists  and  has  a  proper  rational 
transfer  function  matrix  witij  no  poles  iji  Re(s)iO. 
Suppose  that  R,  G(I  +  C  G)  ,  and  S  are  L£e 
stable  and  nonanticipative.  Then, 

(  -G)1  outside  L2e-  Cone  (C,  R,  S)  4Q) 

if  and  only  if  anv  one  of  the  following  conditions 
hold  for  all  real  j) 

il  'min^I+C(-H  C(-jJi))T  ST(  -jJJ)  S<»  (lTC(yi)(G(jX^“ 
-GT(-j'4  RT(-W  R(j*>  *  0  (A4U) 

U)  imin((c(-j*)+c‘(W)T ST(-j*) S(j4(qj‘")+G'1(jn)) 


-RT(-j4  R(j^  2  o 


(A41b) 


iii)  amin(s‘. .w(c(W  +  R  \  j -a)  2  1  (A41c) 

iv)  amax{ R( j4  G(j4^I  +  C(jji)  G(J4)'1  S  (j'Dtj  s  1  . 

'  (A41d) 

When  G*1(jfju)  and  exist,  conditions  (A41a)- 

(A41d)  are  equivalent. 


It  is  trivial  to  show  that  (A41a)  is  always  im¬ 
plied  bv  (A41b)  -  (A4 Id)  and  that,  when  G*  (jui) 
and  R  (jui)  exist,  (A41a)  -  (A41d)  are  equivalent. 

Suppose  that  (A41a)  holds.  Let  (x,y)  by  any 
input -output  pair  satisfying  x  =  -Gy  ;  let 


v  *  S(y-Cx) 


and  let 


v(t)  ,  if  0  *t  <r 
0  ,  otherwise  . 


Let  V'  (ju|  denote  the  Fourier  transform  of  . 
Note  that  from  (A41a),  it  follows  that  for  all 

jViJJi)  I2  -  ,|Rfju>  G(ja> (if  QJ«»  G(jw))'1  S*l(j4 

Vr(j4  !|2  *  0  . 


ilRyl^2  =  ||  RGti+CG)'1^1  v,|T2 
*  II  RG(X+C  G)'1  S'1  vt||t2 
by  the  nonantic  ipativeness  of 

1  a,  G  (U  c  G  f  ,  s~  } 

*JT  1  (  dt 

» 

=  j  :l  R(j*t  G<j4  (i  +  C(j4  GO*))* 1 S"  ij4 


■2 

VJjJ*  ;|  dou 


by  Parseval‘3  Theorem 
and  the  hypotheses  thaj 
(  R.  GII  +  CG)  ,  and  S*  ) 
are  L^e  stable. 
m 

J  ||VT(ju.)H2d-i) 


!|vt',|.2  =  i|v||T2 


=  llS(y-Cx)!lT  •  <A45> 

Conversely,  suppose  that  (A40)  holds.  Let 
Y  and  4)  be  arbitrary.  The  letting  y^Y-e^O* 
ar?d  t-*»°  we  have  from  (A40)  and  Parseval's 
Theorem  that 

l|R(j^0)G(^0)Y0  2  S||S(j»0)(I  +  C(i^)G(j«y)Y0|2 

and  hence  (A41a)  holds.  t  .—I 

This  completes  the  proof  of  Theorem  1. 
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